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Abstract 

We study glueball G production in gluonic penguin decay B —>• G + X s , using the next-to-leading 
order b —>• sg* gluonic penguin interaction and effective couplings of a glueball to two perturbative 
gluons. Subsequent decays of a scalar glueball are described by using techniques of effective chiral 
Lagrangian to incorporate the interaction between a glueball and pseudoscalar mesons. Mixing 
effects between the pure glueball with other mesons are considered. Identifying the /o(1710) to be 
a scalar glueball, we find that both the top and charm penguin to be important and obtain a sizable 
branching ratio for B —> /q(1710) + X s of order 1.3 x 10 _4 (//0.07GeV -1 ) 2 , where the effective 
coupling strength / is estimated to be 0.07 GeV~ 4 using experimental data for the branching 
ratio of /o(1710) —>• KK based on chiral Lagrangian estimate. An alternative perturbative QCD 
based estimation of / is a factor of 20 larger, which would imply a much enhanced branching 
ratio. Glueball production from this rare senri-inclusive B decay can be probed at the LHCb and 
Belle II to narrow down the allowed parameter space. Similar branching ratio is expected for the 
pseudoscalar glueball. We also briefly comment on the case of vector and tensor glueballs. 
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I. INTRODUCTION 


Despite the fact that a glueball state has not been confirmed experimentally, its existence 
is an unavoidable and yet uncannily prediction of QCD. Glueball states may have different 
Lorentz structures in general, such as a scalar, a pseudoscalar or even a tensor with either 
normal or exotic J PC assignment. The prediction for the glueball masses is, however, a 
difficult task. Theoretical calculations indicate that the lowest lying glueball state is a 
scalar 0 ++ state with a mass in the range of 1.6 to 2 GeV. Previous quenched lattice QCD 


Q. 


calculation gave such a glueball mass me equals 1.71 ± 0.05 ± 0.08 GeV |1|. Recent result 
from unquenched lattice QCD calculation |2] gives a result of 1.795±0.060 GeV for this state. 
These results support that the scalar meson /o(1710) to be a glueball. Phenomenologically, 
/o(1710) could be an impure glueball since it may be contaminated by possible mixings with 
the colorless quark-antiquark states that have total isospin zero {3 9]. These mixing effects 
can be either small [s-8] or large js-5, 9], depend largely on the mixing schemes one chose 
to do the fits and these may complicate the analysis. 

Quenched lattice calculation also showed that the lowest lying 0 _+ state for a pure pseu¬ 


doscalar glueball mw 


sum-rule approach 
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rave a higher mass around 2.560 ± 0.035 ± 0.120 GeV [1]. QCD 
131 ] also predicted higher than 1.8 GeV for pseudoscalar glueball 


mass. All these results do not favored the earlier speculation of the 7/(1405) being a pseu¬ 
doscalar glueball. 7/(1405) is indeed a perfect candidate for a 0 _+ glueball since it is copiously 
produced from the radiative decay of J/ip and not seen in the 77 mode. Recently it has 
been shown that when mixing effects and related data for the transition matrix elements of 
anomalous axial-vector current between vacuum and the three states \r]), \r/') and |t/ (1405)) 
are taken into account, the physical pseudoscalar glueball can be as light as 1.4 ± 0.1 GeV 
fid] ] . Clearly, a full unquenched lattice calculation with the fermion determinant included 
must be performed in order to settle this issue satisfactory. Indeed, it has been argued 
some time ago 0 that due to the dynamical fermion effects the full QCD prediction for the 
pseudoscalar glueball mass will be substantially departure from its quenched approximation. 
However, the mass for the lowest lying 0 h state is still missing in the latest unquenched 
lattice results { 2 !]. 

Understanding glueball dynamics is therefore important though very challenging. In this 
work we will study scalar or pseudoscalar glueball production from rare semi-inclusive B rne- 
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son decay. The main purpose of this paper is to refocus the attention of experimentalists for 
undertaking the proposed searches. From the experimental findings we can learn a lot. The 
quantitative numerical predictions are bound to have significant theoretical uncertainties; 
nevertheless experimental efforts will be useful. 


For a recent review for the phenomenology of scalar and pseudoscalar g 
to Ref. 15|, and for a recent summary of glueballs on the lattice, see Ref. 


ueballs, we refer 

Q. 


We lay out this paper as the following. In section II, we discuss the effective vertex of 
b —> sg* induced by the gluonic penguin and the effective couplings between scalar and 
pseudoscalar glueballs with the gluons. In section III, we discuss the interaction of scalar 
glueball with the light pseudoscalar meson octet using the chiral Lagrangian technique. The 
pseudoscalar glueball case will be briefly mentioned as well. In section IV, the rates of 
the semi-inclusive B meson decay into scalar and pseudoscalar glueballs are presented. In 
section V, mixing effects for both scalar and pseudoscalar glueballs will be considered. We 
summarize our results in section VI. 


II. EFFECTIVE COUPLINGS 


Since the leading Fock space of a glueball G is made up of two gluons, production of 


glueball is therefore most efficient at a gluonic rich environment like J/-0 or T —)■ (gg )7 —>■ G 7 
17 . 18j. Direct glueball production is also possible at the e + e _ |l9j and hadron colliders. 
In this work, we show that B —> G + X s decay also provides an interesting mechanism 
to produce and to detect a glueball. The leading contribution for this process is shown in 
Fig. m, where the squared vertex refers to the gluonic penguin interaction and the round 
vertex stands for an effective coupling between a glueball and the gluons. The gluonic 
penguin b —> sg* has been studied extensively in the literature and was used in the context 
for inclusive decay b —> sgr/' 20]. The effective interaction for b —>■ sg* with next-to-leading 
QCD correction can be written as 2lj] 


r^a = ~^^V*V tb s(p') [AF\(g 2 7/i -q,i)L- im b F 2 a^R] T a b(p), (1) 

where AF\ = g)+Ce(q , g))/a s (g) and F 2 = — 2Cs(/i) with Ci(q, g) (i = 4, 6 , and 8 ) 

the Wilson’s coefficients of the corresponding operators in the A B = 1 effective weak Hamil¬ 
tonian, q = p — p' = k + k', and T a is the generator for the color group. We will use the 
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next-to-leading order numerical values of AF± and F 2 2jJ. The top quark contribution gives 
AFi° p = —4.86 and F^° v = +0.288 at /1 = 5 GeV; whereas the charm quark contribution 
involves a q 2 dependence through C^ harrn (q, fi) = Gg harm (g,/i) = P s charm (g,/i) with 

3 charm/„ ^s(aO ~ / >. (10 


and 


\q, +) = (y + Q(q, m c , /i)) , 


Q(q,m, /i) = 4 / dxx (1 — x) In 
Jo 


m 2 — x(l — x)q 2 

H 2 


2 [3 q 2 In(m 2 //i 2 ) — 12m 2 — 5g 2 ] 


9 q 2 

4(2 m 2 + g 2 )\/4m 2 — q 2 

W? 


arctan^ 


4m 2 — q 2 


Here m c = 1.4 GeV is the charm quark mass and C* 2 (+ = 5 GeV) = 1.150. 


( 2 ) 


( 3 ) 



FIG. 1. The leading diagram for B —>• Glueball + X s . 


The following effective coupling between a scalar glueball and two gluons can be 


parametrized as 


22 


£ 0 = fQG^G^, 


( 4 ) 


where Q is the interpolating field for the glueball G, G a ILV is the gluon field strength, and / 
is an unknown coupling constant. When combined with the interaction vertex in Eq. ([T]), a 
scalar glueball can be produced in rare b quark decay, b —> sgG (Fig. (PP>). At the hadronic 
level, this leads to the rare semi-inclusive decay B —* G + X s . The above effective coupling 
has also been used for the study of exclusive decays of B ± —>• / 0 (1710) + 25|. 
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Analogously one can also study inclusive B decays into a pseudoscalar glueball G, whose 


leading effective coupling to the gluon can be parameterized as 


22 


4 , = with G%, =, 


(5) 


where Q denoting the interpolating held for the pseudoscalar glueball and / is an unknown 
coupling. When combined with the interaction vertex in Eq. (JTJ), a pseudoscalar glueball 
can be produced in rare b quark decay, b —$■ sgG (Fig. ([I])). At the hadronic level, this leads 
to the rare semi-inclusive decay B —> G + X s . 

III. DECAY RATES 


If glueballs are composed of only gluons, one can easily obtain the branching ratios for 
B —> G + X s . With the two effective couplings given in Eqs. (JD) and (J3J), the following decay 
rate for b(p) —> s(p')g{k)G{k') (Fig. (JT])) can be obtained readily 


r 


b—>sgG 


'N 2 - l\ G 2 F m 5 b \V*V tb \ 2 ( g- 


4 N c J 

r(i-W ) 2 


x 


dx 


2 5 7T 3 
rv+ 


{m b ff 


47r 2 

dy {|AA’i| 2 c 0 + Re(AFiF 2 *)— + \F 2 \ 2 ^ 


>y- 


y 




with 


y± = 


(1 — x + r) ± \J {1 — x + r ) 2 


4r 


( 6 ) 


(7) 


In Eq. ([ 6 ]), N c is the number of color and 00 , 1,2 are given by 


c o = 2 \~ 2x2 y + (! - y)(y - r )( 2x + y - r )}, 

Cl = (l -y)(y~r ) 2 , 

c 2 = ^[2x 2 y 2 - (1 - y)(y - r)(2xy - y + r)] , 


( 8 ) 


with x — (p' + k) 2 /m 2 , y = (k + k') 2 /m 2 , and r = m^/m 2 . 

The decay rate for pseudoscalar glueball case B —$■ G + X s can be deduced from the scalar 
glueball one by replacing the coupling / with / and the mass me with the pseudoscalar 
glueball mass m& in Eq. (J6]) . This also reproduces previous result for the x — y distribution 
obtained in Ref. 20] for a similar process b —y sgg'. Thus it is not possible to differentiate the 
produced glueball is a scalar or pseudoscalar from the x — y distribution of this semi-inclusive 
B decay. 
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IV. INTERACTION OF GLUEBALLS WITH LIGHT MESONS 


The production rate of a scalar or a pseudoscalar glueball from B decays depend on 
the strength of the couplings / or / respectively. We first discuss how to use the chiral 
Lagrangian methods to incorporate the interaction for a scalar glueball with the other light 
mesons and then employ experimental data to determine /. 

The form of effective coupling of a scalar glueball with two gluons, £o = fQG a lw G a,llJ , 
suggests that glueball couples to the QCD trace anomaly. As briefly mentioned in [26| . 
the interaction of a scalar glueball with light hadrons through the trace anomaly can be 
formulated systematically by using techniques of chiral Lagrangian. The kinetic energy and 


the symmetry breaking mass terms for the light pseudoscalar mesons are given by 

C x = ^[Tr(0*SV„V) + )], 


27] 


(9) 


where f n = 132 MeV being the pion decay constant, £ 2 = £ = exp(2in//„.) with II the 
SU( 3) pseudoscalar meson octet, 

K+ \ 


( J-TT-0 _|_Lit 

y/2* + V6 V 


n = 


7r 


V 


7T 

K~ 


■>° + > K ° 


K 


-feu 


and 


( 10 ) 


X = 2B 0 diag(m, rh, m s ) = diag(m^., m^, 2 m 2 K — m 2 ) 


( 11 ) 


with Bo = 2031 MeV. Here we have neglected the isospin breaking effects due to small 
mass difference between the light u and d quarks and used rh = m u = m^. The QCD trace 


anomaly is well known and given by 


27] 




ba s 

8n 


G^G^ + -£m,qq, 


( 12 ) 


where b = 11 — 2n//3 is the QCD one-loop beta function with rif — 3 being the number of 
light quarks. Treating the effective interaction (j3J) as a perturbation to the energy momentum 
stress tensor, one would then modify to be 

87T 


ba s 

8tt 


(i + f^-e) + E K + f,6) m , 


(13) 
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where f q = fa s m q (16TT\/2/3/3) In [(1 + j3)/(l — f3)\ being the one-loop induced Gqq coupling 
24j with (3 = (1 — 4 m 2 q /rn 2 G ) 1 / 2 . Note that f q is proportional to /. The corresponding chiral 
Lagrangian is thus modified accordingly 


£ * = ^< 1 + /sr c ' Tr 


d^d„T) 


-Jl Tr [e f (x + 2B 0 f x G )e f + e(x + 2B 0 f x g)Z 


(14) 


where f x = diag (f u , f d ,f s ). 

Using the above chiral Lagrangian, one can calculate the decay rates for G —y , 7 r° 7 r°, 


K + K , K°K°, and rfrf. At present there are some data on the ratios of decay width, 


r(7T7r)/r(A'/l) ~ 0.41^0^7 and T{rm)/T{KK) ~ 0.48 ± 0.15 


Theoretically, we obtain 


r(nn)/T(KK) ~ 0.36 and F(r)r))/r(KK) ~ 0.30 which agree with data at about la level. 
With the above chiral Lagragian Eq. ((IT]) , one can also obtain information about the size of 
the coupling / if combined with the branching ratio Br(/o(1710) —>■ KK) = 0.36 ± 0.12 
and the total width T/ 0 (i 7 io) = 137 ± 8 MeV for / 0 (1710) (28]. Using a value of the strong 
coupling constant a s = 0.35 extracted from the experimental data of r decay, we estimate 
the central value of the unknown coupling / = 0.07 GeV -1 . 

In our estimation of / given above, we have extrapolated low energy theorems to the 
glueball mass scale which is higher than Aqcd- The use of chiral Lagrangian is quite ques¬ 
tionable. One might overestimate the hadronic matrix elements in due course. This implies 
that the extracted value of / would be too small. 

An alternative approach is to use perturbative QCD. Indeed, this calculation was carried 
out some time ago for the decay G —>■ KK in 26]. This approach can also give an estimate 
of the amplitude. The problem facing this approach is that the energy scale may not be high 
enough to have the perturbative QCD contribution to dominate. Using the asymptotic light- 
cone wave functions, we find that the branching ratios for G —> KK and G —y m r decays are 
proportional to and /^, respectively. This leads to the ratio of T(ttk)/T(KK) to be about 
0.48 in agreement with data within error bar. Fitting branching ratio for G —> KK, we find 
that the value of / would be about 20 times larger than that obtained above using the chiral 
approach. Incidentally, the value of / derived from the chiral Lagrangian is within a factor 
of 2 compared with that estimated just by using the free quark decay rate of G —> ss. 

In our later calculations, we will use the conservative value / = 0.07 GeV ^ 1 determined 
using the chiral Lagrangian given in Eq. (fT4|) . This will give the most conservative estimate 
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for the branching ratio since the chiral approach gives the smallest /. Clearly, a rigourous 
determination of / is on the lattice. But we do not anticipate to have a lattice evaluation 
of its value any time soon. We have to contend with its rough estimation from all the 
theoretical tools that are available so far. Despite the fact one can not have a precise deter¬ 
mination of the effective coupling /, we believe that glueball production from the proposed 
rare semi-inclusive B decay is a very interesting mechanism which is worth searching for 
experimentally. 

Using the value of / determined above from the chiral Lagrangian, we find the branch¬ 
ing ratio for b —$■ sgG to be 4.5 x 10~ 5 (/ • GeV/0.07) 2 with just the leading top pen¬ 
guin contribution to AF\ is taken into account. The correction from the charm pen¬ 
guin is nevertheless substantial and should not be neglected. Inclusion of both top and 
charm penguins gives rise to an enhancement about a factor of 3 in the branching ratio 
Br(& —> sgG) ~ 1.3 x 10 -4 (/ • GeV/0.07) 2 . Since / 0 (1710) has a large branching ratio into 
KK, the signal of scalar glueball can be identified by looking at the secondary KK invariant 
mass. The recoil mass spectrum of X s can also be used to extract information. The distri¬ 
bution of dBv(b —> sgG)/dMx s as a function of the recoil mass of X s is plotted in Fig. (J2j) 
where M‘j Cg = (p' + k) 2 = mix. 


Recoil Mass Distribution of b -» sgG 



FIG. 2. dBi(B — >• G + X s )/dMx a in unit of 10 5 . In the figure the b quark mass is taken to be 
4.248 GeV, a s = 0.21 at the b quark scale, tb = 1.674 x 10~ 12 s and / = 0.07 GeV -1 . 


Due to parity conservation, a pseudoscalar glueball will not decay into two pseudoscalar 
mesons. Instead it will decay into one scalar and one pseudoscalar mesons or three pseu- 




doscalar mesons. If 7/(1405) is indeed the lowest lying pseudoscalar glueball, kinematics 
allows it to decay into /o(980)t/, KKn, rjTni and so on, as these modes were indeed reported 


as seen by Particle Data Group 28]. The interaction of pseudoscalar glueball with the 


light scalar and pseudoscalar meson octets has been studied in Ref. 30] using the chiral 
Lagrangian techniques for the resonance X(2370) identified as a pseudoscalar glueball. One 
could apply the chiral Lagrangian given in Ref. ^] to the lighter state 7/(1405). However 
the available phase space for 7/(1405) to decay is much smaller than X(2370). We will not 
pursuit that further here. 


V. MIXING EFFECTS 


As is well known, glueballs can mix with other pure qq states with the same quantum 
numbers, which had complicated the identification of glueballs experimentally. For these 
mixing effects, we will adopt the mixing models of Ref. for the^scalar glueballs and 
Ref. [14J for the pseudoscalar glueballs. For the scalar case, we have 


( 


|/o(1370)> 


|/o(1500)> 

V |/o(1710)> J V 
Inverting the above relation, we obtain 


( 0.78 0.51 -0.36 
-0.54 0.84 0.03 
0.32 0.18 0.93 


\ 


( | N)\ 

\S) 

W G )) 


(15) 


\G) = — 0.36|/ 0 (1370)> + 0.03|/o(1500)) + 0.93|/ 0 (1710)) . 


(16) 


Thus the above mixing induces B —* f o (1710)X s , B —> / 0 (1370)X S and B —> / 0 (1500)X S 
with the following relative branching ratios 

Br(/ 0 (1710)) : Br(/ 0 (1370)) : Br(/ 0 (1500)) = (0.93) 2 : (-0.36) 2 : (0.03) 2 = 1 : 0.15 : 0.001 , 

(17) 


up to the corrections from the mass differences among the three states. 
Similarly, for the pseudoscalar case, we have jl^] 


\v) 

W) 

V I?7(1405)) J 


\ ( 0.749 -0.657 0.085 \ / |t/ 8 ) \ 


0.600 0.728 0.331 

-0.279 -0.197 0.940 


l*7i) 

\\G)) 


(18) 
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with 


| G) = 0.085|t 7> + 0.331|r/) + 0.94|?/(1405)) . (19) 

Hence one induces B —y ?7(1405)X S , B —> r/X s and B —> rjX s with following relative 
branching ratios 


Br(i;(1405)) : Br(ij') : Br(ij) = (0.94) 2 : (0.331) 2 : (0.085) 2 = 1 : 0.12 : 0.008 , 

( 20 ) 


up to the corrections from the mass differences among the three states. 

We note that the above relative branching ratios were derived under the assumption 
that B —> (G/G'j + X s arises purely from the three-body process b —>• sg (G/G'j with the 
glueball mixing effects taken into account. There may be two-body process B —>• ( G/G^j + 
X s via b —> s (G/G'j from four-quark operators in the effective Hamiltonian for the G/G 
production from B decay. These contributions will modify the above predictions of the 
relative branching ratios. Being a two-body decay type, the energy of G/G is energetic, 
given by E 2 ~^ dy = (m B + rn 2 G j G — m 2 Xs )/2m B with m B = 5279.50 ± 0.30 MeV and m Xs 
typically equals 1 GeV or so. On the other hand, for the G/G production induced by 
b —>■ sg (G /(5), the energy of G/G is more spread out. One therefore expects with a suitable 
cut on the glueball energy the two-body contributions can be removed, leaving the softer 
G/G production presumably coming from the mechanism we suggested in this paper to 
satisfy the relative branching ratios given in this section. 

To implement the cut, recall that the energy of the glueball E g ^q is related to the rescaled 
variable x — (// + k) 2 /ml according to x = 1 + r — 2 E G / G /mb where r = m 2 G , G /m 2 . With 
E™//, = mb( 1 + r)/2, x min = 0. Thus to make a cut on the maximum of E G / G corresponds 
to a cut on the minimum (lower limit in the integration) of x. In the two-body decay, the 
energy of the glueball is fixed at E 2 ~^ dy = m^l + r — f)/2 with r = m 2 x Jml. Since m Xa 
is of order 1 GeV, we will impose a cut on the minimum of x as in the form x mm — 8 ■ f 
by varying 8 in a range from 1 to 3. The effects of these cuts on the branching ratios of 
Br(T> —>■ /o(1710)X s ) and Br (B —> 77(1405)W S ) are shown in Table [0 The impact of these 
cuts reduce the branching ratios by at most a factor of 3. 
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<5 

a^min — (5 • r 

E G/& = n 2 L { 1 + r ~ xmin ) ( GeV ) 

Br(B^/o(1710)V s ) x 10 5 

Br(B ->■ t?( 1405)X s ) x 10 s 




Top Only 

Top + Charm 

Top Only 

Top + Charm 

0 

0 

2.468 / 2.356 

4.5 

12.9 

7.1 

19.4 

1 

0.055 

2.350 / 2.238 

3.6 

10.2 

5.9 

16.0 

1.2 

0.066 

2.326 / 2.215 

3.4 

9.6 

5.7 

15.2 

1.4 

0.077 

2.303 / 2.191 

3.2 

9.0 

5.4 

14.5 

1.6 

0.088 

2.279 / 2.168 

3.0 

8.4 

5.1 

13.7 

1.8 

0.099 

2.256 / 2.144 

2.8 

7.8 

4.9 

13.0 

2 

0.111 

2.232 / 2.121 

2.6 

7.2 

4.6 

12.2 

2.5 

0.138 

2.174 / 2.062 

2.1 

5.8 

4.0 

10.4 

3 

0.166 

2.115 / 2.003 

1.6 

4.5 

3.3 

8.6 


TABLE I. Branching ratios Br (B —> /o(1710)V s ) and Br (B —> 7(1405)V S ) as function of the 
maximum energy (or minimum x ) cut of the glueballs. We take the nominal value of / = / = 0.07 
GeV" 1 . 

VI. DISCUSSIONS 


Before closing, we would like to make the following comments: 


(1) In 2005, BESII experiment 31] has observed an enhanced decay in J/ijj —>■ 77 j'mr 
with a peak around the invariant mass of 7/7777 at 1835 MeV. Proposal has been 


made to interpret t 
BESIII experiment 


ris state X(1835) to be due to a pseudoscalar glueball [32j. Later 


331 ] has further confirmed the observation of JA(1835) and its spin is 


consistent with expectation for a pseudoscalar. Taking X(1835) to be a pseudoscalar 
glueball, we would obtain a branching ratio of 3.7 x 10" 5 (/ • GeV/0.07) 2 with top 
penguin contribution only, and is enhanced to 1.1 x 1CU 4 (/ • GeV/0.07) 2 if charm 
penguin is also included. If the coupling / is of the same order of magnitude as /, the 
branching ratio for B —> X(1835) + X s is also sizable. 


(2) Lattice QCD calculations also predicted lowest lying spin 2 tensor states: The 2 ++ and 
2~ + masses are 2390 ± 30 ± 120 MeV and 3040 ± 40 ± 150 MeV respectively from the 

n 

earlier quenched calculations , and 2620 ± 50 MeV and 3460 ± 320 MeV respectively 
from recent unquenched results Q]. These heavier tensor states can also be produced 
from the semi-inclusive B decay but with much smaller available kinematic phase space 
compared with the scalar and pseudoscalar states. Leading operators describing the 
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couplings between these tensor glueballs with the gluons are 


g^G; a GT and g^G^fir, ( 21 ) 

where Q ,Ll/ and G^ u are interpolating fields for the 2 ++ and 2 _+ glueball states respec¬ 
tively. Note that Q^ is not the dual of despite the ~ It is straightforward to 
extend the present analysis to these tensor states. 

(3) Lattice calculations [2J] also predicted lowest lying spin 1 vector glueball states which 
would necessarily decay into three gluons. We can also use effective operators de¬ 
scribing interactions between the vector glueballs and gluons. We will leave them as 
exercises for our readers. 

To conclude, we have studied semi-inclusive production of a scalar/pseudoscalar glueball 
in rare B decay through the gluonic penguin and effective glueball-gluon interactions. The 
branching ratio for the scalar glueball case is found to be of the order 10~ 4 . B decays into 
a pseudoscalar glueball through gluonic penguin is also expected to be sizable. Note that 
for the scalar glueball case, we have used a conservative estimate of the effective coupling 
/. Observation of B —» X s / 0 (1710) at a branching ratio of order 10~ 4 or larger will provide 
an strong indication that /o(1710) is mainly a scalar glueball. 

Recall that we have now more than 600 millions of BB accumulated at Belle and more 
than 300 millions at BaBar. Enhancement by a factor of 50 is expected at the Super B 
Factory with a designed luminosity of 8 x 10 35 cm -2 s^ 1 j^. Moreover, at the 14 TeV 
LHCb, the bb cross section can be as large as 500 /ib. One would expect about 50 x 10 12 bb 
pair with an integrated luminosity of 100 fb -1 at the LHCb. Probing the branching ratio 
Br(6 —>■ sgG/G ) at the level of 1CT 4 is quite feasible. One should emphasis that both chiral 
Lagrangian and perturbative QCD approaches used to extract the value of the effective 
coupling / is not precise. Despite this shortcoming, the glueball production mechanism 
suggested in this work is an interesting one. Measurement of this type may also help to 
pin down the energy scale where the chiral Lagrangian is valid for hadronic matrix element 
calculations. We strongly urge our experimental colleagues to rekindle their interests in 
search for the various glueball states; in particular those produced by the semi-inclusive 
decay of the B meson proposed in this work. 
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